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METHOD FOR RECOGNIZING CONNECTABLE SURFACES 



[0001] The invention relates to a method for automatically detecting connectable surfaces in a 
technical system. The system comprises bodies. These bodies can be connected to one another in 
pairs by applying a joining technology. 

[0002] An important joining technology is the production of bonded joints. These are 
increasingly being applied in automobile construction, for example, because it is technically 
impossible to produce a welded connection, the surfaces to be connected are difficult to access 
by the welding operator or automatic welding machine, or because the welded connection cannot 
withstand the loads and forces occurring. A welded connection is often impossible or 
uneconomic in particular, whenever the two bodies are produced from different materials, for 
example, aluminum and steel or aluminum and magnesium, or when at least one of the boundary 
surfaces of the body consists of plastic. 

[0003] The term "joint" also includes below seals between two bodies that have, for example, 
the task of ensuring a minimum spacing between two bodies, and in this case exhibit specific 
elastic properties or effect noise damping or insulation. 

[0004] A computerized design model of the system is prescribed that, for each body of the 
system, comprises at least one surface belonging to the body. 

[0005] The connectable surfaces and the layers between the connectable surfaces are preferably 
decomposed into finite elements. Finite element simulations are subsequently carried out. The 
mechanical behavior of the system is predicted by evaluating the simulation results. 

[0006] The method of finite elements is known from "Dubbel - Taschenbuch fur den 
Maschinenbau", [Dubbel-Manual of mechanical engineering], 20th edition, Springer-Verlag, 
2001, C 48 to C 50, and also from T.R. Chandrupalta & A.D. Belegundu: "Introduction to Finite 
Elements in Engineering", Prentice-Hall, 1991. Strength problems of all types, for example 
relating to stress distribution or stability, are solved numerically by simulation with the aid of 
finite elements. For example, it is determined how a system composed of a number of solid 
bodies is deformed and bent under external loads, and how the bodies are displaced relative to 
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one another. A computerized design model of a system to be investigated is given. A specific set 
of points that are called nodal points are fixed in this design model. The surface or volume 
elements that are formed with the aid of the nodal points as their corners are denoted as finite 
elements. Curved surfaces or bodies that are treated approximately as surfaces, for example 
panels of a motor vehicle bodywork, are often decomposed thereby into shell elements. The 
nodal points form a network in the design model, for which reason the process of fixing nodal 
points and generating finite elements is termed meshing of the design model. Depending on the 
object set, displacements of these nodal points and/or rotations of the finite elements at these 
nodal points or the stresses in these finite elements are introduced as unknowns. Equations are 
set up that approximately describe the displacements, rotations or stresses inside a finite element. 
Further equations result from dependencies between various finite elements, for example because 
the principle of virtual work at the nodal points must be followed, and the calculated 
displacements must be continuous and must satisfy the boundary condition that gaps or instances 
of penetration do not occur in reality. 

[0007] In many cases, such equations are linear in the unknowns. The method of finite elements 
can, however, likewise be applied in the case of nonlinear equations, for example for equations 
in the form of polynomials. Altogether, an often very extensive system of equations with the 
nodal point displacements, nodal point rotations, element stresses or other variables as unknowns 
is set up and solved numerically. The solution describes, for example the deformation state of the 
system under prescribed loads. Stress distributions, vibration behavior, buckling behavior or 
prediction of lifetime, for example, can be derived from this mechanical solution. If, for 
example, displacements and rotations of all the nodal points of a finite element are determined, 
the stress in the element can thus be derived. 

[0008] Various bodies of a system are often meshed independently of one another. For example, 
the system is part of the bodywork of a motor vehicle to be designed, and the bodies are part 
systems that are designed in time-parallel fashion by various suppliers without the meshings 
being adapted to one another. Because the bodies are meshed independently of one another, the 
nodal points on mutually adjoining surfaces of the bodies often do not lie on one another but are, 
for example, displaced relative to one another, or belong to finite elements of different sizes and 
different orientations in space. Such meshings of mutually adjoining bodies are denoted as 
incompatible meshings. 
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[0009] A realistic finite element simulation must take account of the interactions and functional 
dependencies between the bodies that are caused by the mutually adjoining surfaces. What is 
desired are finite element simulations that take account of these interactions and functional 
dependencies even in the case of independent, and therefore generally incompatible, meshings of 
the bodies. The point is that if a compatible meshing were necessary for setting up the system of 
equations and carrying out the simulations, the bodies cannot be meshed independently of one 
another. 

[0010] US 5,560,570 Bl discloses a method for meshing a system with a number of bodies in 
accordance with the finite element method. A computerized design model is respectively 
prescribed for each body. Each body comprises at least one surface belonging to the body. The 
two bodies are firstly meshed independently of one another. The meshing of one body is changed 
in such a way that it matches the meshing of the other body. 

[0011] US 6,343,385 describes a method for determining a trajectory on which a rigid body is 
introduced into a cavity without collision. Both the rigid body and the cavity are meshed in 
accordance with the finite element method. 

[0012] G. Tokar: "Punktschweiflkleber - Eigenschaften und Berechnungsmethode ftir lineare 
Karosseriesteifigkeiten", ["Spot weld adhesives - properties and calculation methods for linear 
rigidity of the bodywork"] VDI- Berichte [Reports] No. 1559, pages 549 - 575, 2000 discloses a 
method for finite element simulation of a bonded joint. Finite element simulations are carried out 
for a system that comprises two sheets that are connected by a bonded seam. Because of external 
loads, displacements occur between and inside the sheets that are predicted by the simulations. 
Finite elements are generated for the simulations in the sheets and in the connecting bonding 
layer. 

[0013] The method disclosed in G. Tokar, loc. cit, requires a great deal of manual work in the 
case when the system to be investigated includes many bodies with connectable surfaces or 
surfaces with complicated geometry. This is the case, for example, whenever the body is a motor 
vehicle to be designed. An operative must mark the bonded seams manually in a design model of 
the system. 
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[0014] It is the object of the invention to provide a method by means of which it is detected 
automatically whether a number of surfaces of various bodies of a technical system can be 
connected to one another by a prescribed joining technology. A computerized design model 
system is prescribed that, for each body of the system, comprises at least one surface belonging 
to the body. i 

[0015] The object is achieved by means of a method as claimed in claim 1. Advantageous 
refinements are specified in the subclaims. 

[0016] In accordance with the method according to the invention as claimed in claim 1, the 
computerized design model of the system comprises a number of surfaces. Each of these surfaces 
belongs to a body of the system. For example, the surfaces are surfaces of the bodies or surfaces 
that approximate the respective bodies. In the case of thin sheets as bodies, the approximating 
surfaces are preferably their middle surfaces. The design model does not necessarily comprise 
volumetric models of the bodies. 

[0017] A joining technology, for example, a specific bonding method, is prescribed. The joining 
technology produces a layer between in each case two bodies, for example a bonded seam or a 
seal. Finite elements are generated for the surfaces. According to the invention, those surfaces or 
sub-areas of surfaces of the system that can be connected by the prescribed joining technology 
are automatically detected. For this purpose, those interspaces between in each case two surfaces 
of the design model that can be filled with a layer produced by the joining technology are 
automatically detected. For example, those interspaces between two surfaces each that can be 
filled by a bonded seam upon application of the bonding method are detected. 

[0018] All the surface pairs that consist in each case of two different surfaces of the design 
model are determined during the detection of the interspaces. Subsequently, connectable pairs of 
finite elements in these surface pairs are automatically determined. The method steps described 
below are carried out for each surface pair, which consists of two surfaces of different bodies. 
The surfaces of such a surface pair are candidates for being connected fully or in sub-areas with 
the aid of a prescribed joining technology. All the element pairs of a surface pair having the 
following properties are selected: 
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- The element pair consists of in each case one finite element of one surface, and of 
one finite element of the other surface of the surface pair. 

- The two finite elements of the element pair have a spacing from one another that 
is smaller than or equal to a prescribed upper bound. 

[0019] An element pair that consists of two finite elements of the same surface is not selected. 
An element pair that consists of two finite elements whose spacing from one another is greater 
than the prescribed bound is likewise not selected. If, for example, the system comprises three 
bodies and the design model of each body in each case comprises four surfaces, and if each of 
these surfaces is decomposed into 100 finite elements, there are thus 4*3/2 = 6 surface pairs 
and 100*100 element pairs per surface pair. If each finite element of one surface has a spacing of 
smaller than or equal to the upper bound relative to four finite elements of the other surface, 
100*4 element pairs are selected per surface pair. 

[0020] This selection is carried out such that all the pairs of connectable finite elements are 
located among the selected pairs, that is to say all the pairs not selected are not connectable. A 
computerized selection rule that can be carried out quickly is applied to the selection. The 
selected pairs of finite elements are thoroughly investigated. It is thereby decided for each 
selected element pair whether the two finite elements of the element pair can be connected by the 
joining technology or not. To take the decision automatically, a decision criterion that can be 
evaluated by computer is applied to compare the positions and/or orientations of the two finite 
elements with prescribed upper and/or lower bounds. These bounds are preferably provided as a 
function of technical properties of the joining technology. For example, in the case of the 
bonding method a bonded seam may be at most 1 mm thick and must be at least 0.2 mm. 

[0021] The invention takes account, without additional method steps, of the possibility that only 
parts of two surfaces can be connected to one another by the joining technology, whereas other 
parts cannot be. For example, one body is a flat sheet and another body is a V-shaped folded 
sheet. One area of the flat sheet can be connected to one limb of the folded sheet, but not to the 
other one. The two sheets are approximated by their middle planes. Connectable finite elements 
of surfaces are determined in accordance with the method according to the invention. In this 
case, it is exclusively finite elements in the one connectable limb of the V-shaped sheet that are 
determined. 



5 



WO2004/109413 



PCT/EP2004/005777 



[0022] Because the tests for connectability are applied to finite elements of surfaces, there are 
fewer comparative operations to be carried out than if the tests were to be applied to finite 
elements in bodies. Specifically, finite elements in surfaces are generally described by fewer 
parameters. The advantage of managing with fewer comparative operations is chiefly important 
when thousands or even hundreds of thousands of finite elements are generated for the surfaces, 
something which can be the case, for example, with design models having many surfaces or in 
the case of a fine decomposition of the surfaces into many small finite elements. 

[0023] A computerized design model of the system with surfaces for the bodies is prescribed for 
carrying out the method according to the invention. There is no need for the design model to 
include volumetric models of the bodies. Consequently, the method can already be applied at an 
early stage in the product development process, specifically at a point in time when only the 
boundary surfaces or approximating surfaces of the bodies, but as yet no details of the bodies, are 
fixed. The use of surfaces and of finite elements in the form of surface elements also effects a 
considerable saving in computing time as well as in arithmetic capability and memory capacity 
by comparison with the use of volumetric models and volume elements as finite elements. 

[0024] Because the pairs of connectable finite elements and therefore connectable surfaces are 
determined automatically, it is impossible for those errors to occur that an operative, for example 
a computing engineer, can commit when manually fixing connectable surfaces. Precisely in the 
case of an extensive system, for example a motor vehicle, many pairs of surfaces come into 
consideration for being connected by the prescribed joining technology. The manual stipulation 
of the pairs that can actually be connected is a time-consuming and error-prone routine operation 
and, in some cases cannot be executed at all in an acceptable time. 

[0025] The decision criterion that is applied in accordance with claim 1 to determine connectable 
element pairs is a computerized, automatically evaluable criterion. It supplies the connectable 
surfaces or areas of surfaces substantially more quickly than an operative by manual stipulation. 
Consequently, the determination of connectable surfaces is objective and practicable and can be 
repeated as often as desired. There is no need to keep asking experienced designers or computing 
engineers for expert knowledge with each application. Subjective factors and errors or mistakes 
that frequently occur in the case of manual stipulation are excluded. There is, furthermore, no 
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need to prescribe a stipulation as to which surfaces are to be neighboring or overlapping. 

[0026] The method according to the invention can also be applied whenever the system 
comprises many bodies with connectable surfaces or surfaces of complicated geometry. For such 
a system it is often impossible to determine interspaces between connectable surfaces by hand in 
an acceptable time. 

[0027] The advantage of automatic detection is still more important whenever the prediction of 
the mechanical behavior must be carried out several times. This is necessary, for example 
whenever various design models of a technical system are to be compared, or various design 
states are run through during designing and, in the process, the positions and/or orientations of 
surfaces are varied. The generation of finite elements is a renewed necessity for each finite 
element simulation of a design model or of a design state. 

[0028] It is possible to determine automatically by means of the method according to the 
invention pairs of connectable surfaces that have not been found by operatives as applications of 
the prescribed joining technology. This is the case whenever finite elements in the surfaces fulfill 
the decision criterion and are determined as being connectable. If, for example, the joining 
technology that is prescribed for the method according to the invention is more cost-effective 
than other joining technologies, the method according to the invention exhibits possibilities for 
savings. For example, bonding is prescribed as joining technology and makes it possible to 
design individual bodies in plastic instead of in steel. It is only by bonding that bodies made from 
plastic can be connected to one another. 

[0029] The method according to the invention can also be applied whenever surfaces of the 
design model have been meshed independently of one another and therefore have incompatible 
meshings. Because the meshings were carried out independently and can be incompatible, the 
bodies of the system can be designed in parallel, for example by different operatives who need 
not synchronize their work for meshing. Because parallel design and parallel meshing is 
rendered possible, and there is no need to coordinate meshings, time is spared and simultaneous 
product design is enabled. It is possible to mesh the surfaces independently of one another and 
firstly to carry out finite element simulations for each body independently of other bodies. Once 
the meshings of the individual surfaces have been produced, they can be reused for a different 
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finite element simulation of the entire system. 

[0030] The selected element pairs detected as connectable delimit interspaces between surfaces 
or sub-areas of surfaces of the design model. Further finite elements are preferably generated for 
these interspaces. In this refinement, the method according to the invention additionally 
facilitates and speeds up the meshing of the connecting layer. The nodal points of these further 
finite elements can be used to set up equations for the mechanical behavior of the layers in the 
interspaces, and for mechanical dependencies between the layers and the adjoining surfaces. 

[0031] The mechanical behavior of a layer can be predicted realistically only when the layer 
occurs as a spatial, that is to say three-dimensional, object, and not as a surface in the simulation. 
Further finite elements are therefore preferably generated for the layer. In accordance with 
claim 12, the interspaces between the pairs of finite elements determined in accordance with the 
invention are automatically meshed. Finite elements with nodal points are thereby generated for 
these interspaces. This meshing need not necessarily depend on the meshing of the 
approximating surfaces. Consequently, the meshing of the layers can be effectively adapted to 
the respective tasks that are to be treated with the aid of the solution of the system of equations 
generated according to the invention. For example, depending on the tasks, the connecting layer 
is decomposed into many small or a few large further finite elements. The thickness of the 
connecting layer is taken into account even when the layer has different thicknesses at various 
points. The layer is treated in the system of equations by using continuum mechanics. For 
example, one body is a flat sheet, and another body is a V-shaped folded sheet. In accordance 
with the method according to the invention it is exclusively finite elements in one connectable 
limb of the V-shaped sheet, and finite elements in the adjoining part of the other sheet that are 
determined as connectable finite elements. Further finite elements are generated only in the 
interspace between the connectable limb and the opposite area of the flat sheet. 

[0032] Mechanical parameters of the layer can also be taken into account in equations of the 
system of equations. If the layer is a bonded seam, for example, mechanical parameters of the 
adhesive used can be taken into account. The mechanical behavior of the connecting layer in the 
event of displacements of the respective surfaces parallel to the layer can be predicted. 

[0033] Claim 2 establishes refinements as to how the selection of element pairs is carried out 
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quickly on the basis of their spacing. The selection that can be executed quickly is made 
according to claim 2 with the aid of the nodal points of the two surfaces of a surface pair. Firstly, 
all the node pairs are determined that consist in each case of one nodal point of one surface and 
one nodal point of the other surface. If one surface comprises N_l nodal points, and the other 
surface comprises N_2 nodal points, NJ * N_2 nodal pairs are determined thereby. The spacing 
between the two nodal points of the node pair is determined for each node pair. A selection is 
made from among the N_l * N_2 pairs of nodal points. Those node pairs whose two nodal points 
have a spacing that is smaller than or equal to a prescribed upper bound are selected. 

[0034] It is possible to determine all the element pairs that respectively consist of a finite 
element of one surface, and a finite element of the other surface, of the surface pair. Very many 
element pairs are often determined thereby. Instead of this claim 2 envisages a preselection: each 
element pair of which one finite element has one nodal point of a selected node pair as a nodal 
point, and whose other finite element has the other nodal point of the pair as a nodal point is 
determined and thereby selected. Further calculations are carried out only for these element pairs 
determined in such a way and thus selected. Those element pairs that have not been selected in 
accordance with the mode of procedure just described are classified as non-connectable. These 
further calculations require calculations that take up more time, as a rule. By contrast, the 
preselection on the basis of the spacings of nodal points can be carried out quickly because 
spacings of nodal points can be calculated quickly. For example, the spacing is determined only 
between selected element pairs, and the element pairs with not too large a spacing are selected 
from among the determined element pairs. 

[0035] Claim 3 and claim 4 develop the refinement according to claim 2. An additional 
preselection from among the determined element pairs is carried out on the basis of the spacings 
of nodal points. 

[0036] In accordance with claim 3, a check is made for each determined element pair as to 
whether each nodal point of one finite element of the element pair has a spacing from at least one 
nodal point of the other finite element that is smaller than or equal to a prescribed upper bound. 
If a nodal point of one finite element has too large a spacing from all the nodal points of the other 
finite element, the test is terminated and the element pair is not preselected and therefore not 
selected and subjected to further tests. Those previously determined element pairs are preselected 
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for which the test supplies a positive result. 

[0037] By contrast, in accordance with claim 4, a test is made for each determined element pair 
as to whether each nodal point of one finite element of the element pair has a spacing from all 
the nodal points of the other finite element that is smaller than or equal to a prescribed upper 
bound. If a nodal point of one finite element has too large a spacing from a nodal point of the 
other finite element, the test is terminated and the element pair is not preselected and therefore 
not selected and subjected to further tests. Those previously determined element pairs are 
preselected for which the test supplies a positive result. 

[0038] Claim 5 provides that the spacing between two finite elements of an element pair is 
compared not only with the upper bound, but also with a prescribed lower one. The element pair 
is not selected if the spacing is smaller than the lower bound. A selection from among the 
element pairs is thereby already carried out on the basis of the spacing. Whenever the spacing is 
greater than an upper bound or smaller than a lower one, it is decided that the finite elements are 
not connectable. 

[0039] Claim 6 establishes refinements as to how the selection of element pairs is carried out 
quickly on the basis of their spacing. Approximations for the spacing are determined thereby 
with the aid of various sequences, and compared with upper and/or lower bounds. At least one of 
these sequences is preferably executed when determining spacing. It is also possible to carry out 
a number of sequences and to compare the respectively determined spacing with an upper and/or 
lower bound in each case. If all the sequences and comparisons lead to a positive result, further 
tests are carried out in order to decide that the two finite elements are connectable. If a 
comparison leads to a negative result at the end of a sequence, it is decided that the two finite 
elements are not connectable. 

[0040] The refinement according to Claim 7 lays down a range of further tests which feature in 
the decision criterion that can be evaluated by computer. At least one of these tests is carried out 
when taking the decision concerning whether the finite elements of a selected element pair are 
connectable or not. The decision criterion preferably applies a logical combination of the results 
of these tests. For example, finite elements of a pair are classified as connectable whenever all 
the tests, or whenever at least a single test are/is fulfilled. The individual tests are preferably 
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carried out in a prescribed sequence such that the individual tests with the lowest computational 
outlay are carried out first. The carrying out of the individual tests is terminated for an element 
pair when it has already been established on the basis of the individual tests already carried out 
whether the finite elements of the pair are connectable or not. 

[0041] At least one of the following individual tests is carried out in accordance with claim 7: 

- Do the finite elements belong to surfaces of different bodies? Specifically, it is 
possible that the two finite elements of an element pair belong to two different 
surfaces of the same body and are connectable. 

The angle between the two finite elements of the element pair is determined, for 
example as an angle between two normals to the finite elements. A test is made as 
to whether the angle is smaller than or equal to an upper bound - the test then 
delivering a positive result - or not. 

One finite element of the element pair is projected along a projection vector. This 
projection vector is generated, for example by generating two normals of the same 
length on the two finite elements, and the projection vector is the sum vector of 
these two (claim 8). It is tested whether the projected finite element overlaps the 
other finite element - the test then supplying a positive result - or not. 

- The midpoints of the two finite elements of the element pair are determined. One 
finite element of the element pair is projected along a projection vector. The 
spacing between the midpoint of the projected finite element and the midpoint of 
the other finite element is determined. A test is made as to whether this spacing is 
smaller than or equal to an upper bound - the test then supplying a positive result - 
or not. 

- As just described, the spacing between the midpoint of the projected finite 
element and the midpoint of the other finite element is determined. The length of 
the longest edge of the two finite elements of the pair is determined. The quotient 
of the spacing and the longest edge length is calculated. A test is made as to 
whether the quotient is smaller than or equal to an upper bound - the test then 
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supplying a positive result - or not. 

[0042] In accordance with claim 9, at least one bound depends on at least one of the following 
parameters: 

- a technical parameter of the prescribed joining technology, 

- the nature of a surface of a body, 

- the material provided for producing a body, 

- a stipulation valid for all the bodies of the system. 

[0043] In the case of a bonded joint, the maximum and the minimum achievable thickness of the 
bonding layer, and the material used for bonding are two such technical parameters. The 
stipulation valid for all the bodies results, for example from esthetic stipulations or from 
company standards. 

[0044] In accordance with claim 10, the term joining technology covers many possible 
technologies, for example bonding, welding or else that of a sealing or insulating or spacing 
layer. For example, a spacing layer made from rubber is inserted in order to observe a prescribed 
minimum spacing between various parts of the bodywork, for example planking and inner parts 
of a motor vehicle. 

[0045] The refinement according to claim 1 1 takes into account the possibility of various joining 
technologies coming into consideration for connecting boundary surfaces. These various joining 
technologies respectively have an evaluation that depends, for example, on the costs and/or the 
reliability of the respective technology. For each pair of boundary surfaces, the joining 
technologies that can be applied for connecting this pair are determined. It is possible that not a 
single, or only one, joining technology is determined. If, by contrast, a number are determined, 
one is selected with the aid of the evaluations. It is possible for different joining technologies to 
be selected thereby for one system. 

[0046] The mechanical behavior of the layer can then be predicted yet more realistically when 
the dependencies and interactions between a layer in one of the interspaces and the surfaces 
connected by the layer are taken into account. Mechanical dependencies exist between nodal 
points of further finite elements of the layer and adjoining points of a boundary surface of a body 
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connected to the layer, for example the principle of virtual work, in accordance with which the 
forces and moments between the nodal points and the adjoining points are in equilibrium. These 
dependencies are taken into account by equations between the nodal points in the layer and 
adjoining points. 

[0047] Claim 1 6 provides an advantageous refinement of how these dependencies are taken into 
account. Claim 17 exhibits a further refinement that spares nodal points and which thereby 
reduces the number of unknowns in the system of equations to be solved. 

[0048] An exemplary embodiment of the invention is described below in more detail with the aid 
of the attached drawings, in which: 

[0049] Figure 1 shows a body and two approximating middle surfaces of the system to be 
investigated; 

[0050] Figure 2 shows surface elements for the body and middle surfaces of figure 1; 

[0051] Figure 3 shows the determination of the minimum and maximum spacing between two 
surface elements (first and second test); 

[0052] Figure 4 shows the fixing of an upper bound for the spacing between two nodal points; 

[0053] Figure 5 shows the determination of the spacing between middle point and point of 
intersection of a normal (first test); 

[0054] Figure 6 shows the determination of the spacing between nodal point and point of 
intersection of a normal (second test); 

[0055] Figure 7 shows the determination of the maximum angle between two surface elements 
(third test); 

[0056] Figure 8 shows the determination of the maximum angle between two surface elements in 
another embodiment (modification of the third test); 
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[0057] Figure 9 shows the determination of the spacing between midpoint and point of 
intersection of a normal (fourth test); 

[0058] Figure 10 shows the determination of the spacing between midpoint and point of 
intersection of a normal and the comparison with an edge length (modification of the fourth test); 

[0059] Figure 1 1 shows the fifth test; 

[0060] Figure 12 shows connectable areas of the surfaces F.l and F.2; 

[0061] Figure 13 shows an interspace that can be connected by a bonding layer; 

[0062] Figure 14 shows an example of two connectable surfaces of the same body. 

[0063] The embodiment described below relates to a bodywork of a motor vehicle as the system. 
The bodywork comprises various panels as well as other bodies, and it is automatically 
determined which of these panels are capable of being connected to one another in which areas 
by bonded joints. 

[0064] A computerized design model of the bodywork was generated with the aid of a tool for 
computerized design (computer-aided design - CAD), and is available in the form of a CAD 
model. The CATIA CAD tool, for example, was used. A description of CATIA is available, for 
example at http://www.catia.com , requested on 2.5.2003. The entire bodywork including the 
panels is designed volumetrically such that the thicknesses of the panels are fixed. 

[0065] Because the panels are very thin by comparison with their extent, they are approximated 
by their middle surfaces in the finite element simulations. All the middle surfaces are 
decomposed in the two-dimensional finite elements in the form of shell elements. 

[0066] A preprocessor is used to generate the data required for a finite element simulation from 
the CAD model for the bodywork. The meshing of the CAD model of the bodywork is carried 
out automatically with the aid of this preprocessor. The method according to the invention which 
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is described below is carried out during the meshing, specifically after the finite elements have 
been generated for the approximating panels. 

[0067] An example of such a preprocessor is the MEDINA software tool. A description of 
MEDINA is available at http://www.c3pdm.com/des/products/medina/documentation/medina- 
DIN4 e.pdf, requested on 2.5.2003. The "MEDINA/PreProcessing" module automatically 
imports a CAD model that is stored in the data format of CATIA or else in the standardized 
STEP or VDA data formats. After the import MEDINA carries out the meshing of the CAD 
model of the bodywork automatically on the basis of stipulations by a user. In this process, the 
finite elements and the nodal points are generated in MEDINA, and these are stored in 
computerized form in the data format of MEDINA. 

[0068] A tool for carrying out a simulation in accordance with the finite element method (FEM 
tool) imports this description in the data format of MEDINA or another data format, and carries 
out the finite element simulations. The person skilled in the art is familiar with various FEM 
tools, for example 

- MSC.NASTRAN and MSC.PATRAN, both described at 
http://www.mscsoftware.com/products/, requested on 2.5.2003, 

- ABAQUS, described at http://www.hks.com/products/products overview.html, 
requested on 2.5.2003, 

- PAMCRASH for finite element simulations of collisions, described at 
http://www.esi-group.com/products/crash/index.php, requested on 2.5.2003. 

[0069] In this example, the design model of the system comprises two sheets and a volumetric 
body K. 1 . The sheets are both 2 mm thick in this example and are approximated in the respective 
middle by two surfaces F.l and F.2. The body K.l is represented by two boundary surfaces F.6 
and F.7. Figure 1 shows the body K.l and four surfaces F.l, F.2 of the two sheets and F.6, F.7 of 
the body K.l. The surfaces F.l, F.6 and F.2 are folded and comprise two planes. The spacings 
and the different orientations of the surfaces and the body in space are represented with great 
exaggeration for the purposes of illustration. The surface F.6 of the body K.l points toward the 
surface F.l and is covered in figure 1. 

[0070] All the pairs of surfaces that belong to two different bodies are determined. In all, there 
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are 4 surfaces and therefore 4*3/2 = 6 pairs, which consist in each case of two surfaces. 
Because the design model comprises a body with two surfaces, one of these six surface pairs 
consists of two surfaces of the same body, specifically the pair (F.6, F.7) The surface pair (F.6, 
F.7) is not investigated for connectability in this embodiment. The remaining five surface pairs 
are investigated. 

[0071] A meshing of all the surfaces is generated. In this example, the finite elements all have 
the shape of triangular or quadrangular surface elements. All four nodal points of a quadrangular 
surface element lie in a plane in this example. Quadrangular surface elements to which this does 
not apply are preferably decomposed into two triangular surface elements for testing for 
connectability. An alternative to this provides that a quadrangular surface element whose nodal 
points do not lie in one plane be replaced for the testing of connectability by an approximating 
quadrangular surface element whose four nodal points all lie in a plane. 

[0072] Figure 2 shows a few surface elements for the bodies and for the middle surfaces of 
figure 1. The quadrangular surface elements preferably have the shape of rectangles, but other 
shapes are also possible. In this example the edge lengths of the surface elements are 10 mm and 
5 mm. 

[0073] The method according to the invention is explained by the example of the two middle 
surfaces F.l and F.6, which approximate two different sheets. The pair of surfaces F.l, F.6 is 
automatically investigated as to which pairs of surface elements of the surfaces F.L, F.6 can be 
connected to one another by one bonded joint each. For this purpose, a decision is taken for each 
element pair as to whether the two surface elements of the element pair can be connected by a 
bonded joint or not. 

[0074] Figure 3 illustrates for example the selection of element pairs for the surface pair that 
consists of the two surfaces F.l and F.6. The nodal points of all the surface elements are firstly 
determined, and their coordinates are stored in one vector each. In the case of five surface pairs, 
therefore, five vectors each having three coordinates of nodal points are stored. The nodal points 
200.1, 200.2, 200.3, 200.4, 200.5 and 200.6 belong to the nodal points of the surface F.l. The 
nodal points 201.1, 201.2, 201.3, 201.4, 201.5 and 201.6 belong to the nodal points of the surface 
F.6. 



16 



WO2004/109413 



PCT/EP2004/005777 



[0075] It is prescribed in this example that a bonded joint may be 1 mm thick at most. An upper 
bound for the maximum spacing between the nodal points of two connectable surface elements is 
derived therefrom. This derivation is illustrated in figure 4. 

[0076] It is assumed by way of simplification in the example of figure 4 that the two surface 
elements 100.2 and 101.2 are parallel to one another. The length of the path from the nodal point 
201.1 of the surface F.6 to the nearest point 230.23 of the surface F.l may be 1 mm at most. The 
nearest point 230.23 is the foot point of a normal to F.6 through the nodal point 201.1. The edge 
lengths of the two surface elements are 5 mm and 10 mm. The spacing 260.1 of the nodal point 
201 . 1 from the nearest nodal point 200.7 is therefore smaller than or equal to 

Va 2 +b 2 +c 2 = Vl 2 + (10 /2) 2 + (5/2) 2 = 5.67 mm. 

[0077] In order to be on the safe side, A_l = 6 mm is fixed as upper bound A_l for the spacing 
between two nodal points. 

[0078] The surface elements of the surface F.l have a total of N_l nodal points, while those of 
the surface F.6 have a total of N_2 nodal points. Each spacing between a nodal point of F.l and a 
nodal point of F.6 is determined. This requires N_l * NJ2 spacing calculations. The calculation 
of the spacing between two points requires much less computing time than other tests of finite 
elements, and so spacing calculations are firstly carried out, and as a function of the result of 
these spacing calculations, element pairs are selected and further tests, requiring more 
computational outlay, are carried out only for the selected element pairs. The calculated spacings 
are buffered in an N_l * N_2 matrix, because each spacing value is used repeatedly. In an 
alternative embodiment, a 1 or a 0 is stored in each of the N_l * N_2 fields of an N_l * N_2 
matrix A. A(i j) is equal to 1 when the spacing between the nodal point No. i of one surface and 
the nodal point No. j of the other surface is smaller than or equal to 6 mm, otherwise it is equal to 
0. 

[0079] In the example of figure 3, the following pairs of nodal points, inter alia, have a mutual 
spacing of at most A_l = 6 mm: 200.1 and 201.1, 200.1 and 201.2, 200.1 and 201.3, 200.1 and 
201.4, 200.1 and 201.5, 200.1 and 201.6. 200.6 and 201.2, 200.5 and 201.3, for example, have a 
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larger spacing. 

[0080] Each element pair is determined whose one finite element has one nodal point of a 
selected node pair as a nodal point, and whose other finite element has the other nodal point of 
the same node pair as a nodal point. A selected node pair in figure 3 is the pair 200.1 and 201.1. 
Consequently, the following 4*4 element pairs are determined whose one finite element has the 
point 200.1, and whose other finite element has the point 201.1 as nodal points: 100.1 and 101.1, 
100.1 and 101.2, 100.1 and 101.3, 100.1 and 101.4, 100.2 and 101.1, 100.4 and 101.1, 100.4 
and 101.2, 100.4 and 101.3, 100.4 and 101.4. 

[0081] A preselection is made from among these determined element pairs on the basis of the 
spacings between nodal points. One of the two following embodiments is applied for this 
purpose: 

[0082] In one embodiment, it is tested for each determined element pair as to whether each nodal 
point of one finite element of the element pair has a spacing from at least one nodal point of the 
other finite element that is smaller than or equal to an upper bound A_2, or not. The bound A__2 is 
fixed such that in the example of figure 4, the element pair (100.2, 101.2) is preselected but the 
element pair (100.1, 101.2) is not. All three finite elements have edge lengths of 10 mm and 
5 mm. As set forth above, however, each nodal point of 101.2 has a spacing from at least one 
nodal point of 100.2 that is smaller than 5.67 mm. Consequently, A_2 = 6 mm is fixed as upper 
bound. 

[0083] The surface element 100.1 has the nodal points 200.1, 200.2, 200.3 and 200.4. The 
surface element 101.2 has the nodal points 201.1, 201.2, 201.3 and 201.4. It is established by 
exercising read access to the N_l * N_2 matrix that the nodal point 200.1 of 100.1 has a spacing 
from the nodal point 201.2 of 101.2 of less than A_2 = 6 mm. Furthermore, it is established that 
the spacings between 200.2 and 201.2, between 200.4 and 201.4 as well as between 200.3 and 
201.3 are also less than A_2 = 6mm. Consequently, the element pair (100.1, 101.2) is 
preselected. By contrast, the spacings between the nodal point 201.3 of 101.2 and the nodal 
points 200.1, 200.4, 200.5 and 200.6 of 100.4 are all greater than A_2 = 6 mm, for which reason 
the element pair (100.4, 101.2) is not preselected. 
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[0084] In the other embodiment, a test is made for each determined element pair as to whether 
each nodal point of one finite element of the element pair has a spacing from each nodal point of 
the other finite element that is smaller than or equal to an upper bound A_2, or not. The bound 
A_2 is fixed such that, in the example of figure 4 the element pair (100.2, 101.2) is preselected, 
but the element pair (100.1, 101.2) is not. All three finite elements have edge lengths of 10 mm 
and 5 mm. The spacing between a nodal point of 100.2 and 101 .2 is at most 

Va 2 + b 2 + c 2 = Vl 2 + 10 2 + 5 2 = 11.22 mm. 

[0085] By contrast, the spacing between 201.1 and 200.3 is more than 12 mm. Consequently, 
A_2 = 12 mm is fixed in this embodiment. 

[0086] In this embodiment, the preselection is undertaken as follows: the surface element 100.1 
has the nodal points 200.1, 200.2, 200.3 and 200.4. The surface element 101.2 has the nodal 
points 201.1, 201.2, 201.3 and 201.4. By exercising read access to the N_l * N_2 matrix, it is 
established that the nodal point 200.1 of 100.1 has a spacing of 10 mm or less in each case from 
201.1, 201.2, 201.3 and 201.4. Furthermore, it is established that the spacing between 200.2 and 

201.1, between 200.2 and 201.2, between 200.2 and 201.3 as well as between 200.2 and 201.4 is 
less than A_2 = 12 mm in each case, that the spacing between 200.3 and 201.1, between 200.3 
and 201.2, between 200.3 and 201.3 as well as between 200.3 and 201.4 is less than 
A_2 = 12 mm in each case and that the spacing between 200.4 and 201.1, between 200.4 and 

201.2, between 200.4 and 201.3, as well as between 200.4 and 201.4 is less than A_2 = 12 mm in 
each case. Consequently, the element pair (100.1, 101.2) is preselected. By contrast, the nodal 
point 201.2 of 101.2 has a spacing of greater than A_2 from the nodal point 200.5 of 100.4. 
Consequently, the element pair (100.4, 101.2) is not preselected. 

[0087] This procedure is carried out for all the determined element pairs. The element pairs are 
selected thereby. The further tests are carried out only for these selected element pairs. 

[0088] The spacing between the two surface elements 100.1 and 101.2 is determined by the first 
test, which figure 5 illustrates. For the test, a normal to the surface element 100.1 is determined 
and a further normal to the surface element 101.2 is determined. A straight line 211.4 is 
generated through the midpoint 240.2 of 100.1. The midpoint 240.2 is determined as the point of 
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intersection of the two diagonals in the surface element 100.1. The straight line 211.4 has the 
same direction as the sum of the two normals to 100.1 and 101.2, respectively. The point of 
intersection 230.1 between the straight line 21 1.4 and the surface element 101.2 is determined. If 
there is no such point of intersection, the test delivers a negative result. Otherwise, the spacing 
between the midpoint 240.2 and the point of intersection 230.1 is compared with a prescribed 
upper bound A_5 = 1 .88 mm. In addition, the spacing is preferably compared with a lower bound 
A_6 = 0.8 mm. If this spacing is smaller than or equal to A_5 and greater than or equal to A_6, 
the test delivers a positive result. Otherwise, it is decided automatically that 100.1 and 101.2 
cannot be connected by a bonded joint. 

[0089] A modification (not illustrated by a figure) of the first test makes provision to determine 
the two midpoints of the two surface elements 100.1 and 101.2. The spacing between the two 
midpoints is determined and used as the spacing between the two surface elements. 

[0090] Figure 6 illustrates a second test for the element pair with the surface elements 100.4 and 
101.1. One normal each is generated to 100.1 at the four nodal points 200.1, 200.4, 200.5 and 
200.6 of the surface element 100.4. The four points of intersection of these four normals with the 
surface F.6 are determined. Such a point of intersection can also lie outside the surface element 
101.1. The normal 210.5 through the nodal point 200.4 and its point of intersection 230.4 with 
the surface F.6 are represented in figure 6. 230.4 lies outside the surface element 101.1. The 
spacing between the nodal point 200.4 and the point of intersection 230.4 of the normals 210.5 
running through 200.3 is determined and compared with an upper bound A_7 and a lower bound 
A_8. Furthermore, 

- the spacing between 200.1 and the point of intersection of the normal running 
through 200.1 with F.6, 

- the spacing between 200.5 and the point of intersection of the normal running 
through 200.5 with F.6, 

- and the spacing between 200.6 and the point of intersection of the normal running 
through 200.6 with F.6 

are determined and compared in each case with A_7 and A_8. Furthermore, four normals that run 
through 201.1, 201.4, 201.5 and 201.6 are generated to 101.1. Their points of intersection with 
F.l are determined. The normal 210.6 through the nodal point 201.1 and its point of intersection 
230.5 with the surface F.l are shown in figure 6. The four spacings of the four points of 
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intersection of the four normals to F.6 with the respective nodal points through 201.1, 201.4, 
201.5 and 201.6 are determined and compared in each case with A_7 and A_8. 

[0091] Figure 7 illustrates the third test, carried out next, by means of which the angle 220.1 
between the two surface elements 100.1 and 101.2 is determined. A normal 210.1 to the surface 
element 100.1, which intersects 100.1 at a foot point 230.10, is generated. In the case of flat 
surface elements, the result of the test is independent of the selection of the foot point 230.10. If 
a surface element with four nodal points is not flat, it is decomposed into two triangular surface 
elements, and the following method is executed for each of these two triangles. The normal is 
preferably of length 1. Furthermore, there is generated to the surface element 101.2 a normal 
210.2 that is likewise of length 1. This normal is displaced to the foot point 230.10. The position 
of this displaced normals is illustrated by the dashed line 210.3. The angle a between 210.1 and 
210.2, which is equal to the angle 220.1 between 210.1 and 210.3, is determined in accordance 
with the following relationship: 

210.1 * 210.2 = | | 210.1 ||*| | 210.2 | | * cos a = cos a = 
cos (220.1) 

[0092] Here, 210.1 * 210.2 denotes the scalar product of the two vectors 210.1 and 210.2, and 
||210. l|| denotes the Euclidean length of the vector 210.1. The angle a determined in such a 

way is compared with a prescribed upper bound A_4 = 10 degrees. The test delivers a positive 
result if the angle 220.1 is smaller than or equal to A_4. Otherwise, it is decided automatically 
that 100.1 and 101.2 cannot be connected by a bonded joint. 

[0093] A modification of the third test just described is illustrated in figure 8. The midpoint 

240.2 of the surface element 100.1 is determined. A normal 210.4 to the surface element 100.1, 
which intersects 100.1 at the midpoint 240.2, is generated. The point of intersection 230.1 of the 
normals 210.4 with the other surface element 101.2 is determined. If there is no such point of 
intersection, the test delivers a negative result and it is decided that 100.1 and 101.2 cannot be 
connected by a bonded joint. If there is a point of intersection 230.1, a normal 210.5 through the 
point of intersection 230.1 is generated on the other surface element 101.2. The angle 220.2 
between 210.4 and 210.5 is compared with the prescribed upper bound A_4. The test delivers a 
positive result if the angle 220.2 is smaller than or equal to A_4. Otherwise it is decided 
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automatically that 100.1 and 101.2 cannot be connected by a bonded joint. 

[0094] The fourth test is illustrated in figure 9. The midpoint 240.2 of 100.1 is determined or 
reused from a prior test. A normal 210.4 to 100.1 is generated through the midpoint 240.2. The 
point of intersection 230.1 of this normals with the surface element 101.2 is determined. The test 
delivers a negative result if there is no such point of intersection. Otherwise, the midpoint 240.1 
of the surface element 101.2 and the spacing between 230.1 and 240.1 are determined. This 
spacing is compared with a prescribed upper bound A_9 = 4 mm. The test delivers a positive 
result if this spacing is smaller than or equal to A_9. Otherwise it is decided automatically that 
100.1 and 101.2 cannot be connected by a bonded joints 

[0095] A modification of this fourth test is illustrated by figure 10. As already described, the 
spacing between the point of intersection 230.1 and the midpoint 240.1 of 101.2 is determined. 
In addition, the length of the longest edge of the two surface elements 100.1 and 101.2 is 
determined. The length of eight edges, specifically if the following edges, is determined for this 
purpose: 

the edge from 200.1 to 200.2, 
the edge from 200. 1 to 200.4, 
the edge from 200.3 to 200.2, 
the edge from 200.3 to 200.4, 
the edge from 201.1 to 201.2, 
the edge from 201.1 to 201.4, 
the edge from 201.3 to 201.2, 
the edge from 201.3 to 201.4. 

[0096] In this case, the edge from 200.1 to 200.4 and the edge from 200.3 to 200.2 are the 
longest edges of 100.1 and are of equal length. The quotient of the spacing between the point of 
intersection 230.1 and the midpoint 240.1 (in the numerator) and the length of the edge from 
200.1 to 200.4 (in the denominator) is calculated. The numerator can vanish, while the 
denominator cannot. This spacing is compared with a prescribed upper bound A_10 = 0.9 mm. 
The test delivers a positive result if this spacing is smaller than or equal to A_10. Otherwise it is 
decided automatically that 100.1 and 101.2 cannot be connected by a bonded joint. 
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[0097] The fifth test is illustrated by figure 11. Two normals 210.1 to the surface element 100.1 
and 210.2 to the surface element 101.2 are formed. The two foot points of the normals can be 
selected as desired. Two normal vectors of equal length are generated on these two normals. 
These two normal vectors are not shown in figure 11. The sum vector 250.1 of these two normal 
vectors is generated. It begins at the foot point 230.10 of the normals 210.1. In this example, a 
straight line 210.8 that runs through the nodal point 200.4 of one surface element 100.1 and has 
the direction of the sum vector 250.1 is firstly generated. This straight line 210.8 intersects the 
surface F.6 at the point 200.13. In the same way a straight line 210.9 that goes through the nodal 
point 200.1 is generated in the direction of 250.1. This straight line 210.9 intersects F.6 at 
200.12. The same procedure is carried out for the two other nodal points of 100.1. A quadrangle 
with the corners 200.12, 200.13, 200.14 and 200.15 is thereby generated. A test is made as to 
whether this quadrangle has an overlap area with the surface element 101.2 or not. If an overlap 
area is present, it is established that the fifth test delivers a positive result. An overlap area is 
present in the example of figure 1 1 . 

[0098] The following tests are preferably carried out for a determined element pair: 

- the modification of the first test (spacing of the midpoints) 

- the third test, if said first test has returned a positive result, 

- the modification of the fourth test, if said third test has returned a positive result, 

- the fifth test, if said fourth test has returned a positive result, 

- if said fifth test has also returned a positive result, it is decided that the two 
surface elements of the element pair are connectable to one another. 

[0099] The following decisions are taken in the example of figure 3 to figure 1 1 : 
100.1 is connectable to 101.2, 

- 100.2 is connectable to 101.3, 

- 100.3 is connectable to 101.4, 

- 100.4 is connectable to 101.1. 

[00100] Figure 12 illustrates which areas of the surfaces shown in figure 1 are connectable to 
one another. These areas are automatically determined by the method described above. It is 
determined which surface elements of F. 1 are connectable to in each case one surface element of 
F.6. The set of these surface elements of F.l delivers the sub-area of F.l connectable to F.6. The 
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corresponding steps are carried out for F.6. One sub-area F.la of the surface F.l and two sub- 
areas F.6a and F.6b of the surface F.6 are shown in figure 12. F.la has the two corner points 
201.15 and 201.16 as well as two further corner points (not shown). F.6a has the four corner 
points 201.11, 201.12, 201.13 and 201.14. The method according to the invention delivers, inter 
alia the result that the two sub-areas F.la and F.6a are connectable to one another. The sub-area 
F.6b is not connectable to a sub-area of F. 1 . 

[00101] The two sub-areas F.la and F.6a of the surfaces F.l and F.6, which are connectable 
to one another are shown in figure 13. The thicknesses of all the sheets of the system are 
prescribed by the computerized design model. Consequently, the two thicknesses d_l and d_3 
are prescribed for those two sheets that are approximated by the surfaces F.l and F.6. Two 
surfaces F.Ik and F.6k are generated. F.Ik lies in the top surface of that sheet which is 
approximated by the surface F.l, and therefore also in the boundary surface of the connecting 
bonding layer. F.6k is congruent to F.6a, but belongs to the bonding layer. F.Ik and F.6k have 
the same dimensions and orientations as F.la and F.6a, respectively. F.Ik lies parallel to F.la 
and F.6k lies parallel to F.6a. The spacing between F.la and F.Ik is 0.5 * d_l (that is to say half 
the thickness of the sheet). Figure 13 shows the two sub-areas F.Ik and F.6k, and also the 
interspace ZW between these two sub-areas. 

[00102] In this embodiment, the two surfaces of a surface pair always belong to two different 
bodies of the system. However, it is also possible to investigate two surfaces of the same body 
for connectability. Figure 14 shows an example for two connectable surfaces F.10 and F.ll of 
the same body. According to the invention, it is determined that a bonded joint can be generated 
that fills up the interspace. 

[00103] In the next step, the interspaces between the connectable sub-areas are preferably 
automatically meshed. The thickness of the sheet is taken into account thereby, and only 
interspaces in layers between sheets are meshed. The meshing of an interspace that connects two 
sheets of the system is executed automatically. The following information is taken over in this 
case from the computerized design model of the system: 

- the spatial position of the two approximating surfaces F. 1 and F.6, and 

- the thicknesses of the two sheets - in this example, each sheet has a thickness that 
is constant over the entire extent, and the two thicknesses can differ from one 



24 



WO2004/109413 



PCT/EP2004/005777 



another. 

[00104] In this example, the thickness of the interspace ZW is 0.8 mm. The thickness and the 
spatial extent of the interspace are obtained automatically from this geometric information about 
the sheets. Instead of this, it is also possible to prescribe the thickness of the interspace and the 
spatial position of the two approximating surfaces. 

[00105] It is possible to decompose the interspaces in the transverse direction into a number 
of volume elements. If, for example, an interspace is 0.8 mm thick and it is provided that an 
interspace in the transverse direction is to be decomposed into two volume elements, volume 
elements are generated that each have an edge length of 0.4 mm in the transverse direction of the 
interspace, that is to say perpendicular to the boundary surfaces of the layer. The meshing of the 
interspaces is controlled by a few parameters, which are clear. These parameters can be selected 
such that the meshing delivers the best results for the respective formulation of the object. The 
volume elements are preferably cuboid, but hexahedrons or other shapes of volume elements are 
also possible. 

[00106] The following prescribed parameters continue to be used for the meshing: 

- a lower and/or upper bound for the edge length of a volume element in each 
longitudinal direction of an interspace, 

- the shape of the volume elements, and 

- a meshing method, for example paving or free meshing. 

[00107] Instead of prescribing an edge length in the longitudinal direction, it is also possible 
to prescribe the number of the volume elements into which the interspace is to be decomposed in 
the transverse direction. 

[00108] All the volume elements preferably have the shape of cuboids or at least of 
hexahedrons. The number of the volume elements in the transverse direction is 2 in this example. 
Thus, in the transverse direction two juxtaposed volume elements are to be generated in each 
case. It is standard for the two volume elements to have the same edge length in the transverse 
direction, and so all the edges in the transverse direction have a length of 0.8 mm : 2 = 0.4 mm. 
Also provided in this example is an edge length in the longitudinal direction of 5 mm in flat areas 
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of the interspace, and 4 mm in curved areas. 

[00109] As an alternative to this, it is not the edge length in the longitudinal direction that is 
prescribed, but a lower and/or upper bound for the ratio of longest to shortest edge of a volume 
element. For example, a ratio of 10 is prescribed in curved areas of an interspace, and a ratio of 
12.5 is prescribed in flat areas. As already explained, the shortest edge length is 0.4 mm. 0.4 mm 
* 12.5 = 5 mm and 0.4 mm * 10 = 4 mm are derived automatically therefrom as length of the 
remaining edges of a volume element in curved and flat areas of the layer, respectively. 

[00110] After conclusion of the meshing of the computerized design model, the physical 
relationships and boundary conditions are supplemented. This step is undertaken, for example 
with the aid of "MEDINA/PostProcessing". 

[00111] An example of such a relationship describes the stress in a finite element as a 
function of the displacement of its nodal points. An expansion tensor s of the finite element is 
determined as a function of the displacement of the nodal points. A compliance matrix D is 
prescribed. The relationship a = D * e exists between the stress tensor a of the finite element and 
the expansion tensor e. 

[00112] It is possible that the deformations result from a temperature variation AT. Let a be the 
expansion coefficient of the material used for producing the respective body. Then the relationship a 
= D * (e - a* AT) exists. 

[00113] Furthermore, the relationship between the acting force F and deformation U is 
determined. A compliance matrix K of the body is derived from properties of the materials that 
are used for producing the respective body, for example elastic modulus and Poisson ratio, and 
from the geometry of the body. The relationship U = K * F exists between the deformation and 
the acting force. It is possible that some components of U are known, for example must vanish, 
and some components of F are known and others unknown. 

[00114] The interspaces are preferably meshed, after determination of the connectable sub- 
areas and of the interspaces between these. However, the meshing is not necessarily executed. 
For example, it is also, possible that instead of this the interspaces are highlighted in the design 
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model. An operative can decide whether precisely these interspaces are actually to become a 
component of a bonded joint, or to be filled up with sealing material and can supplement further 
connectable sub-areas if required, or mark sub-areas detected as connectable as not connectable. 

[00115] It is also possible that the total volume of the interspaces is determined automatically 
and that it is derived therefrom how much material, for example adhesive or sealing material is 
to be filled overall into these interspaces. If a sheet is approximated by one surface, the thickness 
of this sheet is taken into account so that only the volume of the interspace between this sheet is 
taken into account but not the volume of the sheet itself. 

[00116] After the meshing of the boundary surface F.6 of the body K.l, of the middle surface 
F.l of the sheet, and of connecting the bonded joint K.l are concluded, and the system of 
equations has been generated, the system of equations is solved with the aid of a commercial 
software tool for the finite element method (FEM tool). 

[00117] The person skilled in the art is familiar with various FEM tools, for example, 

- MSC.NASTRAN and MSC.PATRAN, both described at 
http://www.mscsoftware.com/products/, requested on 2.5.2003 

- ABAQUS, described at http://www.hks.com/products/products overview.html, 
requested on 2.5.2003, 

- PAMCRASH for finite element simulations of collisions, described at 
http://www.esi-group.com/products/crash/index.php, requested on 2.5.2003. 

[00118] The solution delivers for each nodal point of the design model the value that is 
adopted by the physical quantity at this nodal point. The values of the physical quantity at the 
determined closest points are calculated by substitution in the function. The solution is evaluated 
in order to analyze the design model of the system. 
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[001 19] List of reference symbols 



Symbols 


Meaning 


F.1..F.2 


Approximating surfaces 


F.la 


Connec table sub-area of F. 1 


F.Ik, F.6k 


Boundary surfaces of the connecting layer in the 
interspace ZW 


F.6, F.7 


Two boundary surfaces of the body K.l 


F.6a 


Connectable sub-area of F.6 


F.6b 


Unconnectable sub-area of F.6 


K.l 


Body 


ZW 


Interspace between F.la and F.6a, that is filled up by a 
connecting layer 


100.1, 100.2,... 


Finite elements of the surface F.l 


101.1, 101.2,... 


Finite elements of the surface F.6 


200.1,200.2,... 


Nodal points of finite elements of the surface F.l 


201.1,201.2,... 


Nodal points of finite elements of the surface F.6 


210.1,210.2,... 


Normals to finite elements 


211.4 


Straight line in the direction of the sum vector from two 
normals 


220.1,220.2 


Angle between two normals 


230.1,230.2 


Points of intersection of straight lines with finite elements 


240.1,240.2 


Midpoints of finite elements 


250.1 


Sum vector 


260.1 


Spacing between two nodal points 
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